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Chandra and XMM Observations of the ADC Source 0921-630 

T. R- Kallman 1 , L. Angelini 1 ' 2 , B. Boroson 3 and J. Cottam 1 

ABSTRACT 


We analyze observations of the low mass X-ray binary 2S0921-63 obtained 
with the gratings and CCDs on Chandra and XMM. This object ‘ S 8 
• lirmtinn svstem showing evidence for an accretion disk corona ( ) 

and discuss the implications for accretion disk corona models. For the hrst t 
,nv ADC Lee we use these fits, together with density constraints based 
o”n the O VII line ratio, in order to constrain the flux in the medinm-.on, ration 
Z on rf the ADC Under various assumptions about the source luminosity 
«traLle location of the emitting region. These estimates, together 
with estimates for the emission measure, favor a scenario m w ic 
luminosity of the source is comparable to what we observe. 


1. Introduction 

T ,nite of much study, understanding of the geometry of the gas flows and circumstellar 

. n P x binaries (LMXBs) remains uncertain. The conventional pic ure 

gas in low mass X-ray binar , 1 ' . , arrrption disk is heavily influenced 

consisting of a neutron star, low mass observationa i constraints on 

by analogy with catacl^mt ovmJ* • drive „ s A major contributor 

the accretion disks in LMXBs than mot which makes circumstellar gas a 

mirror for 'x-rajra rather ^han ^medium which can bL studied via its intrinsic emission. In 
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addition, dilution of the circumstellar emission by the strong continuum emanating from 

the neutron star adds to the difficulty of finding signatures of the disk or other structures 
in the spectrum. 

UseM clues to the geometry of LMXBs may be provided by the sources in which 
we are likely observing close to the plane of the accretion disk and binary orbit These 
sources may play a role in our understanding of LMXBs which is analogous to the role 
played by Seyfert 2 galaxies in the understanding of active galax.es; the configuration of 
the circumstellar gas can be studied when the direct X-rays from the continuum source are 
at least partially blocked. In such high inclination systems the presence of partial eclipses 
by he secondary star suggest the presence of an extended component of X-ray emission 
ikely due to scattering and emission from a cloud of highly ionized gas. The X-rays from 
e neu ron star may provide sufficient heating and ionization to account for the properties 

McKiTS'7 dd t 0r0 7h° r n AD< : (Begelman ' Mckee a " d Shields 1982; Begelman and 
p J J II ddlt,on ’ the heated 6^ may flow outward in a wind (Woods et al 1996 

Proga and Kallman 2002) which can potentially affect the mass budget of the system 

The existence of accretion disk coronae has been substantiaUy confirmed, and the column 
density and size have been constrained, by fitting of models to observed X-ray light curves 

al' m i982 m °Ml ‘ToVfil 6 “TIT SPheriCal d ° Ud (White and Holt ' 1982; McClintock, et 

Itone Th7 m u " 7 X ' rayS are refleCted int ° ° ur line of si « ht ^ scattering 
alone These models have also constrained the size and shape of the outer rim of the 

accretion disk, although the results depend on the assumed geometry for the ADC. 

In addition to scattering continuum photons, the corona should also radiate in atomic 
emission features arising from partially ionized material. Such emission has been predicted 

Kalhnan (ImTsT T** (19S9) ’ ^ ** Vr ‘ Uek ' S ° ker a " d Ra y raond (l" 4 ). a"d 
Kal man (1994). Strong line emission is expected under a variety of assumptions, and that 

me emission (as measured by equivalent width) is strongest when the inclination of the 

binary system is closest to 90- Observations of line emission have the potential to provide 

a sensitive test of the hypothesis that the corona is heated and ionized entirely by X-rays 

from the central neutron star, since the ionization balance and emission measure together 

serve to constrain both the gas density and the size of the emission region. 

. Carry , 0Ut SUCh a teSt we have ” sed the Sitings and CCD detectors on 

°n 1 , 6 Two 3 an r "” sate fl' tes to observe the eclipse spectrum of 2S0921-63 a 
well known LMXB which shows the approximately sinusoidal lightcurve indicative of an 
accretion disk corona. Its period of 9.01 days makes it one of the longest period LMXB 
systems, and also makes it particularly suitable for eclipse studies owing to the duration 
of the eclipse. EXOSAT observations of the spectrum through an orbital cycle (Mason 
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et al. 1987) showed that the spectrum softens during eclipse. The spe< ^ u ™ ^ ^ t0 
hard power law, with a photon index of 1.16, and a column density of 1.4x10 (although 
Mason et al. find evidence for contamination at low energies which makes this va ue 
somewhat uncertain). The eclipse (< 50% of the uneclipsed flux) lasted approximately 80 
ksec. The relatively low X-ray to optical flux ratio and many recorded optical dips e.g. 
Branduardi-Raymont et al., 1981, 1983; Chevalier and Ilovaisky, 1981, 1982) are fur er 
evidence for a hidden X-ray source. Cowley et al. (1982) estimate the distance o e ~ 
kpc, implying a luminosity of 2.4 xlO 35 erg s . 

2S0921-63 was observed with ASCA on several occasions, including a total of 27 
ksec of good data obtained during eclipse. Fits to these data reveal a hard power law 
continuum and column density comparable to those observed by EXOSAT. This spectrum 
is significantly harder than the spectra typically observed from low inclination LMXB 
systems. The spectrum contains an iron K line at an energy of 6.75 keV, intermediate 
between the He-like and H-like ionization states of this element. This line is unique among 
LMXBs for its equivalent width; both ASCA and EXOSAT (Asai, 2000; Gottwald et al, 
1996) measure * 120eV for this quantity, while typical LMXBs show much weaker lines 
~ 40 80 eV The flux in the 2-10 keV band in the ASCA spectrum of 2S0921-63 is 3.2 

xio “ erg cm 2 s b In addition, the ASCA spectrum contains evidence for line emission 

in the 0.5-2 keV band, not previously published. 

In the following sections we describe our observations of this source using Chandra and 
XMM, the extraction and analysis procedures, results of model fits, and their implications 
for understanding of 2S0921-63 and for ADC sources in general. 


2. Observations 


2.1. Chandra 

2S0921-630 was observed by Chandra using the HETG and ACIS-S detector between 
UT 15:23:49 2001-08-02 and UT 13:01:13 2001-08-03. The observation lasted approximately 
74000 seconds. In order to best eliminate the influence of any direct X-rays from the 
neutron star, the observation was centered on eclipse phases 0.97 - 1.07, based on the 
ephemeris- A = JD2446249. 18+0.01 + n X9.01 15+0.0005 (e.g. Mason et al., 1987). Although 
the eclipse light curve is gradual, an observation of this duration fits within the time when 
the flux is less than ~ 60% of the uneclipsed value. 

Figure 1 shows the HEG and MEG light curves during our observation, plotted against 
orbital phase using the ephemeris of Mason et al (1987). The data have been binned into 
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100 second bins. A secular variability is apparent, resulting in a factor ~2 increase in the 
counting rate between the beginning and the end of the observation. This variability is 
more pronounced in the MEG than in the HEG, suggesting that the orbital modulation is 
stronger at lower energies. The statistical quality of the HETG spectrum is not sufficient 
to allow sensitive searches for spectral variability associated with this increase. The mean 
first order counting rates were .19 counts/s in the HEG and 0.33 counts/s in the MEG. 
Data was reduced using the standard analysis pipeline and Ciao. We add the + and - first 
order spectra together and analyze only these orders. The zero order ACIS-S spectrum 
showed clear evidence of pileup and will not be discussed further in this paper. The spectra 
from the HEG and MEG were fitted simultaneously, using the xspec package (Arnaud 
1996). Best fit values are found by minimizing the value of the C statistic. Here and in 
the fitting of the XMM data we quote the values of best fit parameters together with the 
errors derived from stepping the parameter in question and determining the A C statistic 
or Ay =10 value. All model fits attempt to account for the ACIS low energy quantum 
efficiency contamination by using the acisabs model in the post-release version of xspec 


2.2. XMM-Newton 

2S092 1-630 was observed by XMM-Newton between UT 11:39:00 2000-12-20 and 
UT 03:24:49 2000-12-21. Data were taken simultaneously for the EPIC detectors MOS 
and PN (using the medium filter) and with the RGS. All data were processed with the 
XMM-Newton Science Analysis Software, SAS 5.3, with the corresponding calibration files 
available for that version. The RGS covers a wavelength range of approximately 6 to 38 A 
(E=0.35 - 2.5 keV) with a resolution of 0.05 A and a peak effective area of about 140 cm 2 
at 15 A (0.826 keV). The observation was separated into two intervals, one lasting 7334 
seconds, and a second lasting 61878 seconds. Here we discuss only the data from the second 
interval. Figure 2 shows the light curve of the PN counts, plotted against orbital phase 
using the ephemeris of Mason et al (1987). The shape of this curve corresponds crudely 
to the shape of the Chandra HETG lightcurve in the local maximum near phase 0.965, 
minimum near 1.01, and overall secular increase. The mean counting rate was .25 counts/s 
in RGS 1 and 0.32 counts/s in RGS 2. The first order spectra from the RGS 1 and 2 are 
fitted simultaneously, using the xspec package (Arnaud 1996). 


3. Results 
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3.1. ASCA 

The results from fitting to the low energy portion of the ASCA spectrum have not been 
previously published, and so we present them here prior to the presentation of the other 
results. The results of fitting the 27 ksec of data contained during eclipse on 04-05-1994 
06:09:12 are shown in Table 1. The spectrum was fit to a cutoff power law plus narrow 
Gaussian emission features. The energies of the features are consistent with those expected 
from the La and L/3 lines of Ne X, and the La lines of Si XIV and S XVI. The O VIII La 
line is not statistically significant in this spectrum. In addition, the fit requires a (sing e) 
iron K line at 1.82 A (6.8 keV), corresponding approximately to the He-like ionization state 
of this element. The measured width of this line is 0.5 keV (Asai et al. 2000), which is 
likely due to the effects of blending. These results foreshadow those found by the other 
instruments described in this paper. 


3.2. Chandra HETG 

In order to better understand the contributions of various components to the spectrum 
we fit the data from all the various instruments with several trial model spectra, with 
progressing physical sophistication. The simplest model consists of a single power law plus 
absorption by neutral material. The results of such a fit to the Chandra HETG data are 
shown in Table 2. Figure 3 shows the ratio of the data to the model for such a fit. T le 
HEG and MEG are plotted separately. The spectra have been binned for the purpose of 
plotting only. The best fit power law index ( 7 ) is 1.01 + oj 5 , which is consistent with that 
found by ASCA. The neutral column is 1.6 xlO 21 cm 2 , which is similar to that found by 
EXOSAT averaged over orbital phase. This value is affected by the treatment of the acis 
contamination in xspec; when the acisabs model is omitted from the analysis the inferred 
column increases to 2.05 xlO 21 cm 2 , other parameters are unchanged. 

Figure 3 shows that, although the fit to the continuum is adequate in the Chandra 
1-15 A (0.8 - 7 keV) region, there are significant residuals near 12, 10.6, 8.4, 7.1, 6.2, 4.7, 
and 1 7-1 9 A (1 1.3 1.44, 1.72, 2, 2.6, 6.4 - 7 keV), which correspond to emission lines 
due to Ne X La, Ne X L/?, Mg XII La and L /?, Si XIV La, S XVI La, and Fe f-XXVk 
respectively. Not all of these lines are statistically significant; Table 2 shows the effect of 
adding multiple narrow Gaussian lines to the absorbed power law model at the energies of 
the significant lines. This results in a reduction of the C statistic from 9390 to 9009 for 
8384 PHA bins All the line fits are consistent with narrow lines (a < 0.005 for the stronger 
lines), corresponding to Doppler widths less than 1500 km/s. The flux in the 2-10 keV band 
in this spectrum is 3.8 xlO 11 erg cm 2 s l . Also given in table 2 and all later tables are 
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the line IDs, and the log of the probability of random occurrence according to an F test. 

The HETG spectrum can also be fitted by a detailed atomic simulation of a 
photoiomzed gas. This was done using models calculated with the XSTAR code (see e g 
Kallman and Bautista, 2000; Kallman et ah, 1996), which assumes that the gas is optically 
thm, and that the ionization is solely due to continuum X-rays from the neutron star. Such 
fits constrain the ionization parameter and normalization to be log(£)~ 4.5 and ~300, 
respectively. Figure 3 and Table 3 show the best-fit model. The statistical quality of the fit 
is significantly better than for the Gaussian fits, C statistic=8731 for 8384 PHA bins. The 
statistical quality of the data is not sufficient to tightly constrain the abundances of the 
emitting ions; formally we obtain 0.1, 0.28, 0.30, 0.33, 0.16, 0.25, 0.22, 2.3e-3, for O, Ne, 
Mg, Si, S, Ar, Ca, Fe, respectively, relative to solar (Grevesse et ah, 1996). 


3.3. XMM 

The XMM data includes the two grating datasets, two EPIC-MOS, and the EPIC-PN. 
Of these, the RGSs provide high resolution (e/Ae ~ 300) over the 0.5 - 2 keV energy band, 
while the EPIC detectors provide moderate resolution (e/Ae ~ 20 - 50) over the 0.5 - 10 
keV energy band. Owing to pileup of the Chandra zero order, the EPIC spectra provide 
our only insight to the broad-band spectrum of 2S0921-63. 


3.3.1. EPIC-MOS 

The MOS provides moderate energy resolution and sensitivity over the entire XMM 
energy band. Figure 5 shows a fit to the global spectrum obtained using a single power 
law with neutral absorption, both the measured and modeled counts and the ratio of data 
to model. Although we have fitted to both MOS detectors simultaneously, only MOS 1 
is shown in this figure for clarity. This shows that the global shape of the spectrum is 
adequately reproduced, and that there are residuals near 6.5-7 keV corresponding to the 
iron K lines, in addition to smaller residuals at lower energies, 0.5 - 2 keV, corresponding 
to the K lines from medium-Z elements. In Table 4 we show the parameters of the power 
law only fit, and for a fit to a power law plus Gaussians at the energies of the iron K 
lines. We defer a discussion of the K lines from medium Z elements to the section on 
RGS fitting. Table 4 shows that the energies of the iron lines are consistent with those 
of the He-like line near 1.82 A (6.8 keV), and the H-like line at 1.78 A (6.97 keV). The 
neutral-like fluorescence line (Fe I-XVII) at 1.94 A (6.4 keV) appears in the residuals but 
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band in this spectrum is 5.9 xlO 11 erg cm s 


3.3.2. EPIC-PN 

The PN provides superior energy resolution and sensitivity at higher energies but 
owing^to calibration uncertainties we can And no satisfactory 

obtained using any of the spectral shapes that we use for any of ^ °ther ^etecto^ 
order to avoid this problem, and still make use of the high energy part of the spectra , 

fi, the PN spectrum only ‘ “ 2 »Ws ft* this Clergy hand. The 

-of the PN spectrum to a simple mode, consisting of 
a power law with index 1.2. This shows residuals at the energies of the 3 iron lines foun 
bv P the MOS and the Chandra HETG, but in addition there is a monotonic decrease in 
flux aWe 7 keV and emission features at energies above 7.1 keV. This spectrum appears 
superficially very similar to the ASCA spectrum of 4U1822-30 in this energy range (Whitd 
1996) This spectrum does not show the behavior expected from Compton reflection w 
is a flattening of the power law at high energies. Instead, the spectrum is appare y 
steepening between 7 and 10 keV and contains narrow features which are statically 
significant. We obtain an acceptable fit this spectrum to a 

a cutoff at 7 keV, Gaussians at the energies of the iron ines. . , • > ^ and 

6 7 and 6 97 keV), plus additional Gaussian emission lines at 1.59 and 1.5 ( • 

« o’kpVl This model fit gives x 2 /i/=874/1028. We have also tested this model against t 
8.2 keV). This mod g X / acceptable. We have tried fitting this 

MOS data in this energy range and found it to be a p ejected 

portion of the spectrum to a single power law plus atomic absorption as m g P 

from an iron K edge or ensemble of edges. The best fit to a single power law con mmun 
plus Gaussian emission at the energies of the iron Kt» lines has absorption e ges a an 
8 4 keV This results in a fit which has x*/^=902/1032, which is significantly greater than 
the at to a curved continuum plus Gaussian emission. The reason for this ,s the presence 
narrow features at 1.59 and 1.51 A (7.8 and 8.2 keV). 

~ r r.zg rsrja rit “ % « 

.in, m m — « »- « f. »v. 

w! can dearly rile out the Ko line of Cu, at f,4 A ,8,5 keV), which is an 
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instrumental feature known to appear in some PN datasets. The results are shown in the 
lower panel of figure 6 and Table 5. The PN data allows a constraint of the widths of the 
ron K lines, for the 1.94 A (6.4 keV) line we can marginally exclude velocities less than 15 
keV, corresponding to Doppler velocities ~ 7000 km s 1, and for the 1.78 A (6 97 keV) line 
we can exclude Doppler velocities less than 430 km s 1 and greater than 3900 km s 1 


3.3.3. RGS 

Chandra HETc'o RGS * similar ^ used for the 

Chandra HETG. One difference is that xspec encounters numerical problems when using 

C statistic with RGS data, so we have rebinned the RGS data using the ftools tool 

grpp a in order to have at least 20 counts per bin, and then used * 2 statistics. The results 

of fitting to various models are shown in Table 6. Figure 7 shows the ratio of the data 

e mo e or the simplest fit, that of a single power law plus neutral absorption The 

spectra have been binned for the purpose of plotting only. The power law index has been 

forced to the value found by the fits to the MOS detectors. For this fit X 2 /» is 1522/1672 

foTh" Yfif 7™ ,Vo eTe r Slgnificant residuals near 22 - 19 > 16 12, 10.6, 10., and 8.6 A 

O VITT T ’ n VTTT T f and 1,44 keV) ’ WhlCh corres P° nd to emission lines due to O VII 1-2 
O VIII La, O VIII 10, Ne X La, Ne X L0, and Mg XII La, respectively. In addition there 

are sigm cant eatuies at 22.9 and 23.1 A which do not have secure identifications. The 
parameters of a fit to the spectrum with Gaussian emission lines added at these energies 
is shown in the second column of Table 6. Inclusion of the narrow Gaussian lines results 
in a reduction in to 1392/1684. Also given in table 6 are the line IDs, and the log 

of p pc ! rand ° m ° CCUrrenCe accord ing to an F test. The spectral resolution 
of the RGS results in upper limits on the widths of the Gaussians which are greater than 

for the Chandra HETG. A notable difference between the Chandra spectrum and that 
shown here is the presence of the He-like O VII lines at 21.8 A and 22.1 A (0.569 and 0.561 
keV). These correspond to the intercombination and forbidden lines, respectively, and their 
best-fit intensity ratio is 1.3. From this we infer that the predominant emission mechanism 
combination, rather than collisional excitation or resonance scattering, since either of 
ese mechanisms is expected to produce strong permitted emission at 21.6 A (0 574 keV) 

In addition, the comparable strengths of the forbidden and intercombination components 

n" K fooTY m the mnge 109 ' 1011 Cm 3 (eg ' Bautista and KalIman 19 "; Porquet and 
aU * greater densities the forbidden component is expected to be collisionallv 

eexci e , w ie at lower densities the intercombination line is weaker owing to the more 
efficient cascades to the upper level of the forbidden line. 
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Table 7 and Figure 8 show the result of fitting the RGS data to an xstar model. We 
aga .in force the continuum slope to match that of the MOS spectrum. We fit to a two 
component model, one component with log(0=l and one with MO=4.5, in order to 
simultaneously fit the low and high ionization lines in the RGS spectrmm The statrs^ 
quality of the fit is significantly better than for the Gaussian fits, X /m- 1418/ 1670. Th 
fit differs from the results of the xstar fits to the Chandra HETG spectra, for which w 
found an acceptable fit to a single ionization parameter, log(0=4.5. An indicator of t e 
difference is the O VII 1-2 line near 21 1 (0.59 keV), since this line requires lower ionization 
parameter in order to be emitted efficiently. Figure 8 shows that this 2 component mode 
adequately accounts for the spectrum in the 15-30 A (0.4 - 0.8 keV) range, but that there 

are inconsistencies between the model and the data in 8 - 15 A (0. • 6 

Most notable is the feature near 9.7 A (1.28 keV) which may be identified with Ne X 

L j3 This feature has a strength which is comparable to the e a 1J | e a ' 

keV) and which is greater than can be produced by recombination or collisions! excitatio . 

ins 11 25 11 45 and 11 7 A (1.148, 1.110, 1.083, and 1.060 
The narrow features near 10.8, 11.25, 11 - 40 , ana n-' ’ . 

keV) are not statistically significant, and so are not fitted by the xstar model. It is also 
worth noting the Ne IX 1-2 lmes near 13.5 - 13.7 A (0.9 - 0.92 keV) which appear to be 
consistent with the predictions of the xstar model. The best-fit abundances for this model 
are: O:Ne:Mg=0.3:L6:2.1 

The differences in the best-fit ionization parameters found from fitting to Chandra 
and XMM can be entirely attributed to the differing bandpasses i of the two > 

The Chandra HETG detects such high ionization lines as S XVI and Si , 

are outside the bandpass of the XMM RGS and which are only emitted a ‘ h '® h 

ionization parameter. The RGS detects the O VII 1-2 lines which are outside the Chandra 
bandpass and which require a lower ionization parameter for efficient emission^ We have 
verified that this is the origin of the discrepancy by fitting both mstruments to only the 
region of wavelength for which both have significant sensitivity, 7- 14 /1 (0.88 - 1.77 key). 
The result is that both fit equally well over a broad range of ionization parameter, since 
both the high ionization and low ionization lines have been excluded in this ht. 


4. Discussion 

4.1. Velocity Offsets and Widths 

In figure 9 we compare the velocity offsets of the lines observed by the various 
instruments The velocity offsets are determined using line wavelengths from the Mb 
database, where available, and from the compilation of Verner (1999) for the others. Figure 
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9 shows that essentially all the line wavelengths are consistent with being emitted at 
zero velocity, although the uncertainties are typically ~1000 km/s or greater. Only the 
wavelength of the Fe XXVI La near 1.78 A (6.97 keV) as measured by the Chandra HETG 

YA?A?DXT tent With Zer ° VeIocity ’ aIthou S h the wavelength of this line as measured by the 
XMM PN is consistent with zero velocity. In the case of the line observed at 1.85 A (6.7 

e ) we include the data for all 3 of the components of the 1-2 lines of Fe XXV, with the 
result that the wavelengths of the forbidden and semi- forbidden lines are consistent with 
smaller velocity offsets than the allows line. This suggests that, if the line-emitting material 

does not have a high bulk velocity relative to us, then the lines are formed predominantly 
by recombination. 

The limits on line widths allow us to approximately constrain the location of the line 
emitting gas, under the assumption that the gas motions are Keplerian in the vicinity of a 
o compact object. If so, the limits on the velocity obtained for the 1.78 A (6.97 keV) 
line using the XMM PN correspond to a radii between 8.8 xlO 8 cm and 7.2 xlO 10 cm. 
Although there are no independent measures of the mass, it is likely that the inclination i is 
close to 90° This can be compared with estimates for the disk outer radius. Using standard 
expressions for the Roche lobe radii (eg. Frank, King and Raine, 1992), and assuming values 
for the mass ratio q in the range 1 - 2.2 (Shabhaz et al., 1999) the radius of the compact 
object Roche lobe is in the range 0.9 - 1.2 xlO 12 cm. Although estimates for the disk radius 
i se are less certain, it is clear that the iron emission originates at a small fraction (~ 10%) 

K In P YA^ bl l 0Uter diSk radiUS ' The Width ° f the 1M A < 6 - 4 keV ) iron Ka line measured 
by the XMM PN is much greater, 0.27 keV. This corresponds to the Keplerian broadening 

a a disk radius of 7 x 10 7 cm. This difference in widths is surprising if it actually reflects 

the range of radn where the lines are emitted. A more likely explanation is that the width 

,° e „ ke ^ kne * s °* ue to blending of multiple components of the fluorescence 

mes from various ionization stages of iron. If so, the estimate of radius obtained from the 
1.78 A (6.97 keV) line remains relevant. 


4.2. Line Luminosities, Variability, and Emission Measures 

The fluxes we derive for the various experiments are 5.4, 3.8, and 5.9 xlO 11 erg 

T ' ’ m the 2 10 k<?V band f ° r ASCA ’ Chandra, and XMM respectively. The fact that 
the ASCA and XMM are consistent, and that the Chandra flux is slightly lower further 
confirms the conclusion that the fact that the Chandra spectrum requires only a high 
ionization parameter component while the XMM RGS requires both high ad low ionization 
parameter components is not due to variability. The lower ionization parameter inferred 
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from the XMM fit would suggest a lower flux during the XMM observat.on if variability 
ter” he origin for the difference. However, the difference is in the oppos.te drrect.on, and 
rr C handra observation has lower flux. From this we conclude that the difference m he 
best Bt ionization parameter for the two observations is due to the differing wavelengt 
bands of the two instruments together with an intrinsic distribution m .on.zat.on parameter 

within the source, rather than variability. 

Figure 11 shows the emissivities for the various lines we observe as a function of 
, c calculated for a library of single-zone xstar photoionization 
IT « we d y o not detect, and whose absence places 

constraints on the conditions in the ADC. These include principally he me of He-1 he 
ions such as Ne IX Mg XII, and Si XIII. From this figure it is clear that the absence ot 
these lines in the Chandra HETG spectra implies ionization parameter values greater than 
f ~ 10* Conversely, the best fit to the XMM RGS spectrum that we find is affected by the 
need to fit the O VII line near 21 A, and this requires ionization parameters ( ~ 10 . 

The line strengths can be used to infer the emission measure of the gas in 2S0921-630 
by dividing the line luminosity by the emissivity of the various lines. As example of such 
calculation is shown in figure 10, for the various strong lines detected by Chandm and X . ^ 
These estimates were derived by using the line emissivit.es shown in figure 11, and summing 
over an arbitrary distribution of ionization parameter binned into 9 bm^ually sf^ 
in log© from 1 to 5. The result is a model luminosity for each line, L, f-M t, 

where , „ is the emissivity for line i in ionization parameter bin k as shown in figure and 

EM k :Z emission melre distribution. This was fit to the observed luminosities and 
errors U‘" and of, derived from tables 2, 5 and 6 using a distance of 7 kpc by varying 
theZnission meakue in each bin independently until a minimum was found in the figure 
of merit We take the figure of merit to be to be L t ’)/' T , ) ■ is proce ure 

is not expected to be highly accurate, owing to the simplified physical 
the coarseness of the various grids, these results in the upper panel of figure 0 show that 
we are able to fit most line strengths to within a factor of 2^ The lower P™ e ot 

shows the derived best-fit emission measure distribution, which corresponds to 3 distmc 

components. This is similar to the result obtained by fitting xstar models to the Chan 
and XMM data separately. This distribution is not consistent with a power law or 

smoothly varying function. 

An equivalent procedure is to use the normalization calculated using the fit to the 
xstar model. The general procedure for doing this is outlined in the xstar manual and in 
the Appendix to this paper. The result is ~ 10” cm which can be compared with 

expected from an ADC. 
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Estimates for the emission measure of an ADC depend on several quantities which 
are not accurately determined, including the geometry of the corona and its effect on the 
ransfer of radiation, the intrinsic luminosity and spectrum of the continuum X-ray source 
the possible role of outflows or winds, and the presence of any other heating or ionization ’ 
mechanism such as magnetohydrodynamic processes (eg. Stone and Miller, 2001) Simple 
estimates can be obtained by assuming the corona is hydrostatic and unaffected by MHD 
processes, as was done by Kallman and White (1989). If so, the corona can be crudely 

represented as a senes of concentric cylinders, each having a Gaussian density distribution 
with height, and where the local scale height is: 


*(*>- 

V GMrriH 


where T is the local gas temperature and M is the mass of the central compact object. The 
maximum density in the corona is attained when the ionization parameter is E = 5* -1-10 

(eg. Krolik McKee and Tarter, 1981), where H = L/(4nR 2 nckT). Then the emission 
measure is approximately 


EM(R) - n 2 max 4-KR 2 z s = 1.2 x 10 61 L 2 s R l0 1/2 T 7 3/2 E 3 

where T 7 is the temperature in units of 10 7 K, R 10 is the radius in units of 10 10 cm, and L 38 
is the continuum luminosity in units of 10 38 erg s 1 . This shows that the emission measures 
inferred from the photoionization model fits can be attained by an accretion disk corona for 
a p ausi e range of parameters, i.e. L 3S - 10 2 , and all other parameters of order unity. 

A further constraint on the ADC comes from the derived density in the O VII lines 
If the intrinsic luminosity of the X-ray source is 10 38 erg s >, then the minimum distance 
at which material at the inferred O VII maximum density of 10 11 cm 3 and ionization 
parameter log(0 ~2 can exist under photoionization is 2.2 x 10 12 cm. This is larger than 
the likely outer radius of the accretion disk, which is less than -10 12 cm (see below). This 
estimate assumes that the local radiation intensity in the ADC is determined by geometrical 
di ution of the X-rays from the compact object. The most likely explanation is that the 
intensity of the ionizing continuum radiation within the O VII region of the ADC is reduced 
by scattering or absorption. If so, the amount of attenuation required in order to reduce 
the O VII radius to a value consistent with the disk is a factor -100. This differs from 
the simple ADC models which have been constructed by eg. Ko and Kallman (1996) and 
imenez- arate et al. (2001), in which ionization parameters low enough to produce O VII 
could only occur at densities - 10 14 cm 3 or greater. An alternative explanation is that the 
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X-ray source is intrinsically faint, ~ 10 36 erg s 1 or less, or that the O VII emission region 
lies outside the traditional estimates for the disk outer radius. 


4.3. Comparison with other LMXBs 

Although 2S0921-630 is not the brightest ADC source, it may the the most suitable 
target for studying the accretion disk corona during eclipse. This is based on the 
spectrum, and is likely due to its long orbital period, which allows long observation 
carried out entirely during eclipse. An illustration of this is provided by the spectrum o 
the ADC source 4U1822-371 (Cottam et al„ 2001), which shows many of the same features 
as 0921-63. However, the line emission from 1822-371 is much weaker in spite of its greater 
total flux In addition, in 4U1822 the emission lines don’t come from the corona. ey 
are phase-dependent and appear to come from the photo-illuminated bulge at the impact 
point This suggests that in 4U1822-371 there may be more dilution of the line emission by 
directly observed X-rays than is expected, or perhaps the inclination is such that we detect 

more direct X-rays than in 2S092 1-630. 


4.4. Neutral-Like Iron line 

The iron K line is unique among X-ray lines from abundant elements in that it can be 
emitted efficiently by neutral or nearly-neutral gas via inner shell fluorescence in addition 
to being emitted efficiently when iron is highly ionized (here and in what follows we 
•neutral-like’ to be ion stages less than Fe XVII, so that the line energy is mdistinguishable 
from neutral). Although many LMXBs (including 2S0921-630) show emission from highly 
ionized iron, a uniformly illuminated accretion disk is expected to radiate a strong line from 
nearly neutral iron. The strength of this line is a constraint on the strength of illumination 
of the accretion disk at radii > 10'cm, where the iron in the disk photosphere is expected 
to neutral-like. A simple estimate for the flux of a neutral-like fluorescence line emitted 

locally from an illuminated disk is: 

J e T h 

where e Th is the threshold energy for photoionization, r FeK is the optical depth at this 
energy and is assumed to be large, Fr Hnuum (R ) '* the continuum flux, w is the fluorescence 
yield and y Fe is the iron abundance relative to solar (Grevesse et al., 1996). we e ne e 

incident continuum flux in terms of the luminosity from the neutron star, L, the radius , 
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a factor / <1 which accounts for non-normal incidence, attenuation between the continuum 
source and the disk, etc., and a factor k which is the fraction of the continuum flux absorbed 
by iron K, we can write: 


L r°° 

Kf(E) MP = 

Then the total fluorescent line luminosity is obtained by integrating over the disk surface 
between two radii and R max : 


J fluorescence — VFe ^ 


L r ft™ 

2 JRmi, 


m 

R 


dR 


If we take f{R)= constant, then 


L fluorescence — Fe^~^ HRm ax /Rmin)f 

Taking plausible values for these quantities: w= 0.34 (valid for neutral-like iron), y Fe = 1, 
k- 0.2 (for an e 1 power law) L = 3 x 10 35 erg s 1 (the measured value for a distance of 
7 kpc), Rmaxj T2 rn ,i rt = 100 (an upper limit with little effect on the result), /= 0.1 (an upper 
limit; cf. Ko and Kallman 1994), we get L fluarescence ~ 5 x 10 33 erg s h This is greater 
than the upper limit allowed by the observations for the 6.4 keV line, 1.8 x 10 33 erg s \ 

suggesting that the illumination is weaker than we have assumed, / < 0.03, or that the iron 
abundance is less than solar. 


5. Summary 

The complementary capabilities of Chandra and XMM provide new insight into the 
reprocessing gas in 2S0921-630. The XMM RGS allows detection of line emission from 
elements such as oxygen, while the Chandra HETG allows detection of the K lines of sulfur 
and iron. The XMM EPIC detectors allow measurement of the spectrum above 7 keV and 
measurements of line widths. We have found line emission from medium-Z elements, O, 
Ne, Me, Si, S, and also K lines from iron. With the exception of the neutral-like and H-like 
components of the Fe K line as measured by the PN detector on XMM, the lines have 
widths which are less than can be measured with the Chandra and XMM detectors. The 
measuied widths for the two components of the iron line suggest upper limits to the size 
of the ADC of 7xlO I0 cm for a 1 M© compact object at 90° inclination if orbital Doppler 
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broadening is the dominant mechanism. The line spectra from each grating mstrumen 
adequately with photoionized models. However, the XMM observations require a second 
component at lower ionization parameter than inferred for Chandra^ The O V M are 
detected in the XMM RGS and imply gas densities in the range 10 - 10 cm , whic 
in turn implies that either the X-ray source is intrinsically faint, that the emission region 
lies outside of the accretion disk, or that large regions of the ADC are very optically m 
The distribution of emission measures of xstar models which comes closest to fitting the 
ensemble of the line data is a 3 component model, similar to that inferred from the separate 
xstar fits to the Chandra and XMM grating data. The emission measures in erre rom 
both a simple optically thin DEM analysis and from the xstar fits are consistent wi 
simple estimates for an ADC if the intrinsic luminosity of the source is approximate y 
we observe. The high energy spectrum as observed by the PN shows compkx structure at 
energies above 5 keV : curvature in the continuum spectrum, multiple iron Ko components 
similar to those seen in the Chandra HETG, plus either additional emission features at 7.8 
keV and 8.2 keV, absorption, or reprocessed emission which differs from standard reflection 

models. 

The picture of the ADC that we are left with is as follows: The total ADC emission 
measure is comparable to that obtained from simple estimates, although the distribution 
of ionization parameter is not smooth, monotonic, or easily predicted or compared with 
models. The width of the Fe XXVI La line implies that high ionization gas is distributed 
primarily at radii less than cs 7 x 10‘°cm. The intensity of the neutraUlke Fe Ka line 
implies that this material is primarily at large radii, i.e. greater than ~ cm. , 
measured width of this line is due to blending of multiple components due to different 
stages. The constraints on the density of the O VII emitting gas require ^at either the 
intrinsic luminosity of the continuum source be low, i.e. < 10 erg s , or that the O 
emitting region be heavily shielded from direct X-rays. The first scenario implies that we 
observe the full luminosity of the continuum source. This conflicts with ADC mode s l in 
which the intrinsic luminosity of ADC source is much greater than we observe and that 
much of what we observe is scattered or difused through the ADC. The second scenario 
requires a large column density (> 10“ cm ’) of absorbing or scattering gas which is not 
apparent from the observed spectrum. The high energy spectrum, above the iron K line 
region shows evidence for complicated and unexplained emission, or possibly absorption. 
Although simple models for absorption do not fit the spectrum as well, both emission and 
absorption models produce fits which are statistically acceptable. The absorption scenario 
is more plausible owing to the likelihood of dense gas associated with the disk atmosphere, 
bulge, or wind impinging on the line of sight. If so, this gas is highly ionized and also has 
large column density (> 10 24 cm 2 ). More sensitive observations, or more secure E 
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calibration, in the 1-20 keV energy band may allow this to be tested; an absorption scenario 
would predict that the underlying power law would recover at energies -15 keV. 
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Appendix 

The general procedure for inferring the emission measure from the normalization calculated 
using the fit to the xstar model is as follows. The observational data is fitted to a library of 
xstar table models constructed using the xstar2xspec tool. When this library was calculated 
a source luminosity L xstar was used, and when the data was fitted to the model a best fit 
normalization, xstar ionization parameter (0 and column density (N) were derived The 
best-fit xstar model was calculated as a spherical shell, and its emission measure is 

EM xstaT = n 2 4wR 2 AR = 4tt— atarjV 

If the distance to the source (D) and its intrinsic luminosity (L) are known, then the 
emission measure of the source is given by 

EM = normalization x EM xstar 
where D^ pc is the distance to the source in kpc. 

In the case of the fit to the Chandra HETG spectrum the best fit normalization is 210 
for an ionization parameter £ ~ 10 5 , and L xstar = 10 6 and EM xstar = 2 x 10 61 . This ’ 
corresponds to a source emission measure of EM = 6.8 x 10 54 cm 3 L J D; where L, 0 is 
the continuum luminosity in units of 10 38 erg s 1 and D 7 is the distance in units of 7 kpc. 
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Fig. 11.- Emissivities in erg cm 3 s 1 for strong lines in the Chandra/XMM band as 
calculated by xstar 
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Table 1. ASCA Lines 


quantity value 


Absorption, NH (xlO 22 cm 2 ) 
Power Law Index (photons) 
Cutoff Energy (keV) 
normalization (erg cm 2 s *) 
Flux 2-11 keV (erg cm 2 s 1 ) 
line 1 

Wavelength (A) 

Energy (keV) 

normalization (xlO 5 erg cm 2 s 
Equivalent Width (eV) 

Line ID 
log(P) 
line 2 

Wavelength (A) 

Energy (keV) 

normalization (xlO 5 erg cm 2 s 
Equivalent Width (eV) 

Line ID 
log(P) 


0.17 

1.1 

17.8 

4.0+0. 1 xlO 3 
3.45 xlO 11 

12.4 

i n+0.02 
0.75 

J ) 9.7+6°' 

22+8 
NeX La 
-4.3 

10.7 

i ic +0.05 
110 0.08 

*) 4.7+f; 

22+8 
Ne X L/J 
-3.0 
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Table 1. ASCA Lines, continued 


quantity 


value 


line 3 

Wavelength (A) 

Energy (keV) 

normalization (xlO 5 erg cm 2 s ) 
Equivalent Width (eV) 

Line ID 
log(P) 
line 4 

Wavelength (A) 

Energy (keV) 

normalization (xlO 5 erg cm s ) 
Equivalent Width (eV) 

Line ID 
log(P) 
line 5 

Wavelength (A) 

Energy (keV) 

normalization (xlO 5 erg cm 2 s ! ) 

Line ID 

log(P) 

Equivalent Width (eV) 

X 2 /PHA bins 


6.2 
1.99+0.1 
3.5+?; 
20+11 
Si XIV La 
- 2.2 


4.77 
2 . 6 + 0. 1 
1.5+2 
15+9 
S XVI La 
-0.32 


1.82 
6 . 8 + 0. 1 
4.7+1;? 
Fe K 
-3.49 
115+40 
718/595 
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Table 2. Results of fitting to Chandra HETG 


Parameter 

NH (xlO 21 cm 2 ) 
photon index 

normalization( erg cm 2 s 1 ) 
line 1 

wavelength (A) 
energy (keV) 

normalization (xlO 3 erg cm 2 s 
equivalent width (eV) 

ID 

log(P) 
line 2 

wavelength (A) 
energy (keV) 

normalization (x 10 5 erg cm 2 s *) 
equivalent width (eV) 

ID 

log(P) 
line 3 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s ! ) 
equivalent width (eV) 

ID 

log(P) 
line 4 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s 
equivalent width (eV) 

ID 

log(P) 

line 5 

wavelength (A) 
energy (keV) 

normalization (x 10 5 erg cm 2 s 
equivalent width (eV) 

ID 


Power Law Only Power Law + Gaussian 


1.6 

1.6 

L 00 + o ° i 5 

1 - 09 + oi 5 

• 0031+;™ 7 

- 0033 +X 

- 

12.16 

- 

i no+0.0015 

0.00025 

— 

2 4+ 2 

1.5 

— 

6.65 

- 

NeX La 

— 

- 19.2 

- 

10.60 

- 

1 . 17 + 0.01 

- 

1 1 + 1-2 
11 0.8 

— 

3.52 

- 

NeX L0 

■ 

- 7.74 

- 

8.56 

- 

1 . 45 + 0.0025 

- 

7-0+5 

- 

0.185 

- 

Mg XII La 


- 0.06 

- 

6.18 

- 

2 . 006+°°°° 3 

- 

1 - 95+0 2 

— 

10.2 

- 

Si XIV La 

- 

- 34.3 


4.73 

2.62+0.01 
1.03+} 5 
13.2 

S XVI La 


UoVP\ 


a no 
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Table 2. Results of fitting to Chandra HETG, continued 


Parameter 


Power Law Only Power Law + Gaussian 


line 6 

wavelength (A) 
energy (keV) 

normalization (x 10 5 erg cm 2 s ) 
equivalent width (eV) 

ID 

log(P) 
line 7 

wavelength (A) 
energy (keV) 

normalization ( x 10 5 erg cm 2 s ) 
equivalent width (eV) 

ID 

log(P) 
line 8 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm s ) 
equivalent width (eV) 

ID 

log(P) 

C statistic/PHA bins 


9390/8384 


1.938 

6.400+0.01 

1 . 1+2 

26.2 

Fe neutral-like Ka 


-2.03 


1.865 

6.650+0.01 

1.4+2 

31.1 

Fe XXV 1-2 
-3.28 


1.797 

6.900+0.01 
< 2 
50.1 

Fe XXVI La 


- 0.2 

9009/8384 
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Table 3. 


Results of fitting to Chandra HETG with xstar table + power law model 


Parameter Value 


NH (x 10 21 cm 2 ) 2.0 

photon index 1.03+og 5 

normalization (erg cm 2 s .00286+ ^ 

xstar column density 10 21 cm 2 

lo g(0 4.5 

xstar normalization 300 . 0 + 2 $) 

C statistic/PHA bins 8731/8384 
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Table 4. Results of fitting to XMM EPIC MOS 


Parameter 


Power Law only power law + Gaussian 


NH (xlO 21 cm 2 ) 
photon index 

normalization (xlO 3 erg cm 2 s 
line 1 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm s ) 
equivalent width (eV) 

ID 

log(P) 
line 2 

wavelength (A) 
energy (keV) 

normalization ( x 10 J erg cm 2 s ) 
equivalent width (eV) 

ID 

log(P) 

X 2 /*' 


1.735+0.05 

1.757+0.05 

1.184+ 0.01 

1.193+ 0.01 

6.7+0.15 

6.78+0.15 


1.85 



6.698+0.1 



4. 0+1.0 



1500 

_ 

Fe XXV 1 - 2 

- 

-7.5 


1.78 



6.978+0.1 



1. 2+1.0 



288 



Fe XXVI La 



-1.95 

2077/1209 

1969/1203 
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Table 5. 

Results of fitting to XMM EPIC PN 

Parameter 

Power Law Only 

Power Law 4- Gaussian 

NH (xlO 21 cm 2 ) 

6.6 

11.2 

photon index 

1.2 

1.2 

High Energy Cutoff (keV) 

_ 

7 

normalization ( erg cm 2 s 

*) .0112+0.005 

.0116+. 0003 

line 1 

wavelength (A) 

— 

1.94 

energy (keV) 

— 

6.38+0.1 

width (keV) 

- 

0.29+0.15 

normalization (xlO 5 erg cm 

2 S 1 ) ~ 

1 4+0.15 

± n fit: 

equivalent width (eV) 

- 

u.uo 

14 

ID 

- 

Fe I - XVII Ka 

log(P) 

— 

-99 

line 2 

wavelength (A) 

— 

1.86 

energy (keV) 

- 

6.68+0.03 

width (keV) 

- 

<0.1 keV 

normalization (xlO 5 erg cm 

2 s>) - 

5. 7+0.2 

equivalent width (eV) 

- 

63 

ID 

- 

Fe XXV 1-2 

log(P) 

- 

-5.59 
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Table 5. Results of fitting to XMM EPIC PN, continued 


Parameter Power Law Only Power Law + 


Gaussian 


line 3 

wavelength (A) 
energy (keV) 
width (keV) 

normalization (xlO 5 erg cm 2 s x ) 
equivalent width (eV) 

ID 

log(P) 
line 4 

wavelength (A) 
energy (keV) 
width (keV) 

normalization (xlO 5 erg cm 2 s x ) 
equivalent width (eV) 

ID 

log(P) 
line 5 

wavelength (A) 
energy (keV) 
width (keV) 

normalization (xlO 5 erg cm 2 s x ) 
equivalent width (eV) 

ID 

log(P) 

x 2 / u 


1.78 

o.9o o.ot 
0.05+0.04 
4.0+0. 1 
49 

Fe XXVI La 
-99. 

1.59 
7.8+0. 1 

< 0.5 
1.3+0. 3 
23 

-3.74 

1.51 

8 . 2 + 0. 1 

< 0.3 
1.4+0. 3 
25 

-5.58 

2074/1039 874/1028 
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Table 6. Results of fitting to EPIC RGS 


Parameter 

NH (xlO 21 cm 2 ) 
photon index 

normalization (erg cm 2 s *) 
line 1 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s ! ) 
equivalent width (eV) 

ID 

log(P) 
line 2 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s 2 ) 
equivalent width (eV) 

ID 

log(P) 
line 3 

wavelength (A) 
energy (keV) 

normalization (x 10 5 erg cm 2 s *) 
equivalent width (eV) 

ID 

log(P) 
line 4 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s M 
equivalent width (eV) 

ID 

log(P) 


Power Law Only Power Law + Gaussians 


2.1 

1.2 

.0085+0.0001 

23.22 

0.534+0.01 

5.3+7 

2.8 

-1.06 

22.92 

0.541+0.01 

7.1+8 

3.5 

-1.34 

22.10 

0.561+0.01 

4.5+5 

2.2 

O VII 1s 2 ( 1 S)-1s2s( 3 S) 
-1.15 

21.77 

0.569+0.01 

5.9+12 

3.3 

O VII ls 2 ( 1 S)-ls2p( 3 P) 
-2.27 


2.1 

1.2 

.0085+0.0001 
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Table 6. Results of fitting to EPIC RGS, continued 


Parameter Power Law Only Power Law + Gaussians 


line 5 

wavelength (+ 
energy (keV) 

normalization (xlO 5 erg cm 2 s *) 
equivalent width (eV) 

ID 

log(P) 
line 6 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s *) 
equivalent width (eV) 

ID 

log(P) 
line 7 

wavelength (A) 
energy (keV) 

normalization (xlO 5 erg cm 2 s x ) 
equivalent width (eV) 

ID 

log(P) 
line 8 

wavelength ( A ) 
energy (keV) 

normalization (x 10 5 erg cm 2 s *) 
equivalent width (eV) 

ID 

log(P) 

line 9 

wavelength (A) 
energy (keV) 

normalization (x 10 5 erg cm 2 s ! ) 
equivalent width (eV) 

ID 

log(P) 

X 2 /" 


1522/1672 


19.0 

0.652+0.01 

7.6+2 

5.4 

O VIII La 
-14.35 

15.97 

U-MO 0.00025 

1.4+2 

1.2 

O VIII Lp 
-0.84 


12.13 

i noi + a0003 

0.004 

3 . 6+11 

3.8 

NeX La 
-2.17 


10.65 

1 1 £7+0.0003 
0.004 

2.0+0; 5 

2.9 

Ne X L/3 
-0.60 


8.40 

1.474+g;882 3 

4.0+I:? 

7.3 

Mg XII La 
-3.09 

1393/1650 
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Table 7. Results of fitting to XMM RGS with xtar table + power law model 


Parameter 

Value 

NH (xlO 21 cm 2 ) 

2.2 

photon index 

1.2 

normalization (erg cm 2 s *) 

•0081 + ;£*£ 

xstar column density 

10 21 cm 2 

xstar density 

10 11 cm 2 

Component 1 
log(0 

1 

xstar normalization 

2.0x10 4 

Component 2 
log(0 

4.5 

xstar normalization 

26.9 

x 2 / u 

1413/1670 








